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A . Abstract 

D : 

^ . Tau-pair production in the process e+e -^ e~^e t^t was studied using data 

collected by the DELPHI experiment at LEP2 during the years 1997 - 2000. 

k>( ; The corresponding integrated luminosity is 650 pb~^. The values of the cross- 

;_i ' section obtained are found to be in agreement with QED predictions. Limits 

^ on the anomalous magnetic and electric dipole moments of the tau lepton are 

deduced. 
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1 Introduction 

This paper presents a study of tau pair production in photon-photon colhsions using 
the data collected by the DELPHI detector at LEP in the period from 1997 to 2000 
(LEP2) at collision energy ^/s between 183 and 208 GeV. The total integrated luminosity 
used in the analysis is 650 pb~^. At LEP this process was first observed by the OPAL 
collaboration [1] and subsequently studied by the L3 collaboration [2]. 

The final state e"'"e~r+r^ can be produced via a set of Feynman diagrams. In this 
paper we present the cross-section measurement for the contribution of the so-called 
multiperipheral graph (Fig. 1) which corresponds to collisions of two virtual photons. 
The same final states produced via other diagrams (less then 1% of the cross-section) are 
considered as a background. 



Figure 1: The dominant diagram for the reaction e"'"e — > e"'"e t^t 

The study of the reaction e"''e~ -^ e^e~T^T^ explores two fundamental problems. First 
of all it provides a deep test of QED at the level of the fourth order in a. Furthermore, 
the 7rr vertex is sensitive to the anomalous electromagnetic couplings of the tau lepton. 
Since the multiperipheral e"'"e~ -^ e"'"e~r^r~ process diagram contains two such vertices, 
the anomalous magnetic and electric dipole moments can be extracted by comparing the 
measured cross-section with QED expectations. 

The rest of the paper gives a detailed description of the e"'"e~ -^ e"'"e~r+r~ cross- 
section measurement, including tau-pair selection, background estimation, selection and 
trigger efficiency calculation and systematic error estimation. In the last part of the paper 
the measured cross-sections are used to derive limits on the anomalous electromagnetic 
moments of the tau lepton. 

2 Monte Carlo simulation 

The signal process was simulated using the Berends, Daverveldt and Kleiss generator 
RADCOR (BDKRC) [3], which calculates the cross-section for the multiperipheral dia- 
gram with radiative corrections on the electron and positron lines. The following signal 
definition was used: the invariant mass of tau pairs had to be less than 40 GeV/c^; both 
beam particles had to be scattered by less than 10 degrees; and at least one of them had 
to be scattered by less than 2 degrees. With these restrictions the accepted cross-section 
was 1.44±0.04% lower than the total cross-section predicted by BDKRC (which is about 
450 pb at LEP2 energies) for the unrestricted phase space. The r decay was simulated by 



the TAUOLA package [4], which includes photon radiation from the decay products. The 
BDKRC generator was also used to estimate the background coming from the process 
e"'"e~ — s> e^e' ji^ji^ . 

To simulate the e"'"e~ -^ e+e~e"'"e~ background, the Berends, Daverveldt and Kleiss 
generator DIAG36 (BDK) [5] was used. Hadron production in two photon collisions was 
simulated by PYTHIA 6.1 [6]. Non-multiperipheral four-fermion processes (such as WW, 
ZZ, Zee and others) were simulated by WPHACT [7]. 

The generated events were passed through the full simulation program of the DELPHI 
detector and were reconstructed with the same program as for the real data [8]. 



3 Event selection 

In most events produced by two-photon collisions both beam particles scatter at small 
angles and remain undetectable inside the beam pipe. Therefore only the decay products 
of the tau leptons can be seen in the detector. To suppress background, only one-prong 
decay channels with one tau decaying into an electron and the other into a non-electron 
(hadron or muon) were considered. The analysis was based entirely on the measured 
tracks of charged products of tau decays; the neutral particles from tau decays were 
ignored in this analysis. 

To select runs with good performance of the sub-detectors [8]- [9], only runs with the 
Time Projection Chamber (TPC), the Forward Chambers (FCA, FCB) and one of the 
additional barrel tracking detectors (ID or VD) fully operational were retained. Table 1 
presents the luminosities used in the analysis, luminosity-weighted centre-of-mass energies 
and energy ranges. 





1997 


1998 


1999 


2000 


Luminosity, pb^^ 
< Ecm >, GeV 
Energy range, GeV 


52.3 

182.7 
182.7 


152.6 
188.7 
188.7 


224.2 

197.6 

195.5-201.5 


217.5 

206.3 

204.5 - 208.0 



Table 1: The integrated luminosities, mean collision energies and collision energy ranges. 



The event selection procedure was divided into two steps. The first step (preselection) 
selected a sample of two-photon events with two good tracks which were not back-to- 
back in the plane perpendicular to the beam axis. A track was considered as good if 
the momentum derived from its curvature was greater than 100 MeV/c, momentum error 
better than 100%, polar angle 9 between 20° and 160°, and impact parameter with respect 
to the interaction point below 10 cm along the 2;-axis ^ and 5 cm in the r — plane. 

The following cuts were applied in this first selection step: 

• There had to be exactly two good tracks from particles with opposite charges, at 
least one of them having momentum greater than 300 MeV/c. 

• To suppress background from fermion pair production, the total energy of two 
charged particles had to be less than 30 GeV. 



^The DELPHI coordinate system has the 2:-axis aUgned along the electron beam direction, the x-axis pointing toward 
the centre of LEP and the j/-axis vertical, r is the radius in the {x, y) plane. The polar angle d is measured with respect 
to the z-axis and the azimuthal angle (/) is about z. 



• To enrich the sample with e+e -^ e+e r^r events, the acoplanarity of two tracks 
^ had to be greater than 0.5° and their resultant transverse momentum greater than 
500 MeV/c. 

• To select events with a high trigger efficiency, the transverse energy, defined by 

Et = El sin 6i + E2 sin 62 , 

where Ei and E2 are the energies of the two charged particles and 61 and 62 are 
their polar angles, had to be greater than 2 GeV. 

• In the year 2000, the operation of one of the twelve TPC sectors was unstable and 
the dE/dx measurement vital for this analysis was poor, so events with at least one 
track in or near (closer than 10° in 0) to this TPC sector in 2000 were rejected. 

• Finally, to ensure the transverse momentum balance of 77 system, single and double 
tagged events were rejected by requiring that no energy deposition in the forward 
electromagnetic calorimeters (STIC or FEMC) exceeded 60% of the beam energy. 

The last cut suppressed the events with highly virtual photons. About 90% of events 
passing this cut had the momentum transfer — g^ less then 1 GeV^/c^. After applying 
the cuts described above, the predicted event composition in the preselected sample was 
as follows (1999 data): 

e+e" -^ e+e-e+e" 41 % 

e+e" -^ e+e-/i+;U- 47 % 

e+e- -^ e+c-r+r- 8 % 

e+e^ — s> e+e^qq 3 % 

e+e" -^ T+T- 1 % 

The fraction of other events was less than 1%. The efficiency of the preselection for 
e"'"e~ — s> e^e^T^T^ events was of the order of 5%, the largest suppression of the signal 
coming from the requirement of exactly two good tracks seen in the detector (about a 
factor of 4) and from the cut on Et (about factor of 2). Figure 2 shows the comparison 
between data and simulation of the distributions of invariant mass, total energy, total 
transverse energy and total transverse momentum of the pair of charged particles. The 
e"'"e~ — s> e^e~T^T~ events are shown by the shaded histogram. The Monte Carlo is 
normalised to the luminosity of the real data. The data deficit is mainly due to the 
trigger inefficiency which is corrected at the later stages of analysis. 

In the final step of the selection, the event was retained if one of the charged particles 
was identified as an electron and the other as a non-electron. This step was based on the 
TPC measurement of the di?/dx pulls for the muon, electron, kaon and proton hypotheses. 
The dE/dx pull for a specific particle hypothesis is defined as the ratio 

{dE/dx)raeas-{dE/dx)exp ,,. 

iix = , Uj 

O'AE/dx 

where {dE /dx)exp is the value expected for the particle X with given momentum and 
o'AE/Ax is the error of the measured energy loss (d^/dx)meas- To check the dE/dx calibra- 
tion, test samples of e^e~ -^ e~^e~fi^fi~ and e+e~ -^ e"'"e~e^e~ events were picked out 
from the preselected sample. A small angular dependence of the dE/dx measurements 
was found as well as some disagreement between data and simulation. Corrections which 
were functions of azimuthal and polar angle were applied to the measured dE/dx values. 
Residual disagreement was removed by scaling and smearing the specific energy loss mea- 
surement in the simulated events. Independent calibrations of real and simulated data 



"J 



Acoplanarity is defined as 180° — 102 — <^il- 



were performed for each year of data taking analysed. The efficiency to measure dE/dx 
is discussed later in the paper. 

With corrected dE/dx information, a track was identified as an electron if 11^ > 3 and 
as a non-electron if lie < —3. Figure 3 illustrates the particle identification cuts. The 
distributions of the pulls for the electron and muon hypotheses are shown for the 1999 
real data and simulation. Each distribution is shown after applying all selection cuts 
except the cut on the variable shown. The shaded histograms show the signal. 

A considerable amount of kaon and proton background from e^e~ -^ e^e~qq events 
remained after the cuts on the pulls for the muon and electron hypotheses. Figure 4 (left) 
shows the specific energy loss for electron candidates plotted versus the momentum of 
the particle. Proton and kaon bands are clearly visible. To remove the kaon and proton 
background, the electron selection was tightened. The dE/dx for the electron candidate 
had to not exceed 1.9 times the minimum ionisation, and the pulls for the proton and 
kaon hypotheses for the electron candidate both had to be outside the ±1.5(T interval: 
IIIkI > 1.5 and lllpl > 1.5. Figure 4 (right) shows the distribution of the pull for the 
proton hypothesis with all selection cuts applied except the cut on the variable shown. 
The hatched histogram shows the background from e^e~ — > e"'"e~qq events, the shaded 
histogram shows the rest of the background. The cuts on this variable are indicated by 
arrows. 

Table 2 summarises the efficiency of the preselection, final step of selection and overall 
selection efficiency. The drop in the preselection in 2000 is caused by the removal of 
events in or near the unstable TPC sector. The uncertainties of the selection efficiency 
determination are discussed later in the paper. 





1997 


1998 


1999 


2000 


preselection 
selection 


5.39 
17.3 


5.37 
16.4 


5.38 
16.4 


3.85 
16.1 


overall 


0.93 


0.88 


0.88 


0.62 



Table 2: Efficiencies (%): preselection, final step of selection and overall efficiency. 

In total 2390 candidate events were selected. Figure 5 compares the distributions 
of electron and non-electron candidate momenta for selected events to the simulation 
prediction for the combined 1997-2000 data. Figure 6 shows the visible invariant mass 
distribution for selected events for the same data sample. Trigger efficiency is taken into 
account in these distributions (see below). 



3.1 Trigger efficiency 

The low momenta of the r decay products in the process e"'"e~ -^ e"'"e~r"'"r~ and the 
presence of only two tracks in the event could make the probability of triggering on such 
an event considerably below 100%. The determination of the trigger efficiency is therefore 
important in this analysis. 

The trigger efficiency was estimated from the subsamples of selected events using 
the fact that an event can be detected by different components of the DELPHI trigger 
system [10]. Trigger subcomponents of the tracking system were combined into barrel 
and end-cap triggers. For events with one track in the barrel and the other in the end-cap, 
the number of events detected by the barrel trigger (Nb), by the end-cap trigger (Ne), 



and by both triggers (Nbe) were counted using the decision functions of the trigger. 
The barrel and end-cap single track trigger efficiencies were calculated, for electrons and 
non-electrons separately, by the formulae: 



N, 



^Barrel 



BE 



iV, 



Nf 



^end- 



BE 



cap 



N, 



(2) 



B 



Finally, the efficiency of the DELPHI calorimetric trigger to the whole event was estimated 
from the events triggered by any of the tracking detectors using a similar technique. The 
results of the trigger efficiency calculation are summarised in Table 3. The track pair 
trigger efficiency was calculated from the "OR" of the single track efficiencies using the 
ratio of the barrel and forward tracks predicted by simulation. The tau pair trigger 
efficiency was calculated as "OR" of the tracking and calorimetric triggers. 





1997 


1998 


1999 


2000 


Barrel track 














electron 


71.4 ±17.1 


94.4 ±5.4 


84.6 ±7.1 


92.3 ±7.4 




non-electron 


100t?7.9 


85.2 ±6.8 


85.0 ±8.0 


78.6 ±11.0 


End-cap track 














electron 


26.3 ±10.1 


36.5 ±6.1 


21.5 ±4.6 


22.0 ±6.9 




non-electron 


31.3 ±11.6 


30.4 ±6.1 


25.9 ±4.8 


23.1 ±5.8 


Track pair 




94.5+?-? 


95.5 ±2.7 


93.3 ±3.5 


93.5 ±4.1 


Calorimetry 




6.7 ±1.9 


8.6 ±1.2 


7.1 ±0.9 


7.7±1.1 


Tau pair 


94.91^+ 


95.9 ±2.5 


93.8 ±3.3 


94.0 ±3.8 



Table 3: Summary of the trigger efficiency measurements (%) 



3.2 Efficiency of the dE/dx measurement 

Both tracks of the selected event had to have specific energy loss measurements. An 
imperfect detector simulation can lead to a discrepancy in the dE/dx measurement effi- 
ciency for good tracks in real and simulated events. To take into account this possible 
disagreement, the efficiency for a good track to have a dE/dx measurement (to be a 
"good TPC track") was calculated for e+e" -^ e"'"e~e"'"e~ and e+e~ -^ e'^e'fi^fi^ sam- 
ples extracted from preselected events (efficiencies to measure dE/dx of pions and muons 
were assumed to be equal). Muon events were selected by requiring at least one track to 
be identified by the muon chambers and electron events were selected using information 
from the DELPHI RICH detectors. For muon and electron events the efficiency to be a 
"good TPC track" was determined from the ratio 



N, 



e 



dE/dx 



dE/dx 



(3) 



tot 



where N^E/dx is the number of events with both tracks having a dE/dx measurement and 
Ntot was the total number of selected events in the given sample. The tau-pair efficiency 
was estimated as the product of the single track efficiencies for muon and electron. The 
tau-pair efficiencies derived from e+e" -^ e'^e~fi~^fi~ and e+e~ -^ e+e~e"'"e~ events for 
data and Monte Carlo are presented in Table 4 and were used for selection efficiency 
correction and for systematic error estimation. The selection efficiency was multiplied 



by the factor '^^^''"'....^^ and half of the correction was included into the systematic error 
together with the uncertainties from the test sample statistics. 





1997 


1998 


1999 


2000 


Efficiency in data, % 
Efficiency in MC, % 


82.9 ±1.0 
82.3 ±0.4 


82.6 ±0.6 

82.5 ±0.2 


82.4 ±0.5 
82.3 ±0.1 


83.5 ±0.6 

84.6 ±0.2 



Table 4: Summary of "good TPC track" efficiency estimations. These efficiencies are 
already included in the total efficiency in Table 2. 



3.3 Residual background 

Several sources of background for e^e" 



e^e t'^t events have been considered: 



• The background from e^e~ — > e^e~qq, mainly from protons and kaons selected due 
to the tails of the dE/dx pulls for the proton and kaon hypotheses; 

• The background from e+e^ — * e"''e"e+e^ and e+e^ — > e~^e~fi'^fi~ events due to the 
tails of the distributions of the dE/dx pulls for the electron and muon hypotheses; 

• Background due to other four-fermion processes: non-multiperipheral diagrams (in- 
cluding e"'"e~r+r~ final states) and multiperipheral process e+e~ -^ e"'"e~r+r~ which 
does not satisfy signal definition; 

• The process e^e~ — ^ t^t^ (background from other fermion pair production processes 
was found to be negligible). 

The background fractions for the main background sources and their uncertainties are 
summarised in Table 5. The contribution from other background sources was negligible. 
The theoretical precision of e+e~ -^ e^e~qq generation by PYTHIA is not well known. 
Therefore it was estimated from the real data by inverting the dE/dx cut on the electron 
candidate: dE/dx >1.9 M.I.P. instead of dE/dx <1.9 M.I.P. After comparing these test 
samples enriched with e+e^ -^ e^e~qq events with the simulation, an error of 20% was 
ascribed to the e"'"e~ -^ e"'"e~qq event generator. 



Channel 


1997 


1998 


1999 


2000 


e"'"e~ — s> e"'"e~qq 


4.3 ±1.5 


3.2 ±0.7 


3.3 ±0.8 


3.2 ±0.8 


e+e" -^ e+e'e+e" 


2.7 ±0.2 


3.4 ±0.1 


4.0 ±0.1 


2.4 ±0.1 


e+e" -^ e^e~fj,^ 


2.9 ±0.1 


5.0 ±0.1 


2.4 ±0.1 


3.8 ±0.1 


other 4-fermion 


1.5 ±0.3 


1.5 ±0.3 


1.2 ±0.2 


1.2 ±0.2 


e+e" — ^ TT 


0.69 ±0.01 


0.55 ±0.01 


0.47 ±0.01 


0.40 ±0.01 


Total 


12.1 ±1.5 


13.6 ±0.8 


11.4 ±0.8 


11.0 ±0.8 



Table 5: Summary of background fractions. The numbers are the expected fractions (%) 
of the specified backgrounds in the selected sample. Errors are statistical errors of the 
simulated samples and theoretical uncertainties of the Monte Carlo generators added in 
quadrature. 



4 Systematic error estimation 

The following sources of systematic error on the cross-section measurement were con- 
sidered: uncertainties of selection and trigger efficiencies and uncertainty of background 
level. Track selection, event selection and the statistical error of the simulated samples 
were taken into account in the calculation of the uncertainty in the selection efficiency. 

The systematic error arising from track selection was estimated in the following way. 
Each cut of the track selection was varied typically by the size of the resolution of the 
corresponding variable from its nominal value in both directions. The corresponding 
change of the cross-section A was compared to the value of the expected statistical 
fluctuation a due to the non-identical event sample. If the value A was less than a, no 
systematic error was ascribed to the corresponding cut; in the opposite case the value 
of VA^ — cr^ was included in the systematic error. The systematic error arising from 
varying the event selection cuts was estimated in a similar way. 

To calculate the systematic error due to the angular corrections applied to the dE/dx 
measurements, the dE/dx correction functions were varied by the uncertainty of each of 
their parameters and the analysis chain was repeated. The variation of the measured 
cross-section was added to the systematic error. The systematic errors corresponding to 
scaling and smearing the pulls were calculated similarly. 

The systematic errors associated with track selection cuts, event selection cuts and 
dE/dx corrections are summarised in Table 6. The numbers are given for 1999 data (for 
other years uncertainties of most of the sources scale approximately as inverse square 
root of the statistics). Additional contributions to the selection efficiency uncertainty 
also presented in Table 6 come from the statistical error of the Monte Carlo sample and 
the selection efficiency correction described in Section 3.2. 

The largest contribution to the systematic error is given by the uncertainty of the 
trigger efficiency determination, dominated by the statistics of the real data events, see 
Section 3.1 and Table 3. An additional contribution arises because the trigger efficiency 
for background events, assumed to be equal to that of the signal, may be different. A 
conservative estimate of this uncertainty was obtained by changing the trigger efficiency 
for background upwards to 100% and downwards by the same amount. 

The systematic error due to residual background includes the simulated sample statis- 
tical uncertainty and the theoretical uncertainty of the Monte Carlo generators, mainly 
for the e+e~ -^ e+e~qq process, see Table 5. 

The sources of selection efficiency uncertainty are described in detail in Table 6. 

The sources of systematic uncertainty are summarised in Table 7. Total systematic 
errors calculated as the sum in quadrature of all described components are also presented 
in Table 7. The following uncertainties were assumed to be fully correlated between differ- 
ent years: generator theoretical error; trigger efficiency for background; and uncertainties 
estimated from variation of track and event selection cuts. Systematic errors from other 
sources were treated as uncorrelated. 



5 Results of the cross-section measurement 

The cross-sections were computed using the formula 

Nobs — Nbg 

^ = ^ ^ r (4) 



syst. error source 




value, % 


track selection cuts 








J^imp 


0.7 




^imp 


1.1 




6p/p 


0.7 


event selection cuts 








He 


0.3 




n^ 


0.3 




acoplanarity 


0.6 


dE/dx corrections 








^eO 


1.0 




ne0 


0.9 




n^^ 


1.0 




n^.^ 


1.0 




scaling 


0.7 




smearing 


0.6 


MC statistics 




0.8 


"Good TPC track" 






correction 




0.6 


Total 




3.0 



Table 6: Systematic errors for 1999 data coming from track selection, event selection, 
dE/dx corrections, simulated samples statistics and "good TPC track" correction. 



where Nobs is the number of observed events, N^g is the expected number of background 
events in the assumption that background events have the same trigger efficiency as signal 
events, Ssd is the selection efficiency, Strig is the trigger efficiency and C is the integrated 
luminosity. 

The numbers of observed and expected events, the measured cross-sections and the 
cross-sections from the BDKRC Monte Carlo simulation together with their ratios are 
presented in Table 8. The predicted number of events was calculated from the signal 
and background simulation, taking into account trigger efficiency and corrections to the 
dE/dx efficiency. Agreement was found between the measurements and the Standard 
Model (SM) predictions calculated by BDKRC. The ratio of observed and predicted cross- 
sections was averaged over all LEP2 data, taking into account correlations of systematic 
errors. The result was found to be 0.96±0.04. The average LEP2 cross-section is 429±17 
pb corresponding to the luminosity- weighted mean centre-of-mass energy of 197.1 GeV. 
The cross-section predicted at this energy by BDKRC is 447.7 ± 0.3 pb. 



6 Determination of anomalous magnetic and electric 
dipole moments 

In the Standard Model, leptons are considered as point-like objects. Therefore the 
observation of a deviation of the magnetic or electric dipole moments of the leptons from 
their SM values would open a window onto the physics beyond the SM. The anomalous 
magnetic moments of the electron [11] and muon [12] are known with high precision. 





1997 


1998 


1999 


2000 


Trigger eff. 
Selection eff. 
Background 
Luminosity 


7.0 
5.1 
1.7 
0.6 


2.7 
3.2 
0.9 
0.6 


3.6 
3.0 
0.9 
0.6 


4.5 
3.0 
0.9 
0.6 


Total 


8.9 


4.3 


4.7 


5.4 



Table 7: Relative systematic errors on cross-section (in 



Year 


Observed 


Expected 


'^measj PU 


(Tmc, pb 


O'meas/O'MC 


1997 


211 


224 ± 18 


401 ± 32 ± 36 


428.2 ±0.5 


0.94 ±0.11 


1998 


629 


652 ± 24 


419 ±19 ±18 


436.7 ±0.5 


0.96 ±0.06 


1999 


909 


937 ± 39 


436 ± 16 ± 21 


448.5 ± 0.5 


0.97 ±0.06 


2000 


641 


665 ± 32 


443 ± 20 ± 24 


459.4 ±0.5 


0.97 ±0.07 



Table 8: The numbers of observed and expected events, measured cross-sections, QED 
predictions and their ratios. The first error on the measured cross-sections is statistical, 
the second is systematic. 



but the short life-time of the tau-lepton does not allow measurement of its anomalous 
moments with similar precision by a spin precession method. 

The generalised form of the tt'-/ vertex can be parametrised as follows: 



-ieu{j>'){Fi{q^)^^ + iF2{q^)a 



^u 



+ Fsiq'h'a 



-}uip)ef,{q) 



(5) 



2mr 2mr 

where e^(g) is the polarization vector of the photon with momentum q. The form factor 
Fi describes the distribution of electric charge and e^- = eFi{0), while F2 and F3 are form 
factors related to the anomalous magnetic moment a,- and electric dipole moment dr'- 

«r = ^V^ = ^2(0) (6) 



and 



i^3(0) 



2mrdr 



(7) 



In the SM at tree level, a,- = and dr = 0. Accounting for loop diagrams gives a non-zero 
value to Qr = 11773(3) ■ 10~^ [13], while a non-zero value of dr is forbidden by both T 
invariance and P invariance. 

The values of a,- and dr have been measured by several groups. The L3 and OPAL 
collaborations [14,15] studied radiative Z -^ tt'J events and set the following 95 % CL 
limits on the values of the anomalous magnetic and electric dipole moments: 



-0.052 <ar< 0.058 
-0.068 <ar < 0.065 



and I^T-I < 3.1 
and I^T-I < 3.7 



1-16 



^10 "" e ■ cmj 
'10^^^ e ■ cm) 



(L3) 
(OPAL). 



The best limit so far on dr was obtained by BELLE [16]: 



-16, 



-0.22 < ?R:e{dr) < 0.45 (10~ -e ■ cmj 
-0.25 < Qm{dr) < 0.08 (lO^^^e ■ cm) 
Other limits on a^- and dr can be found in [17]. 
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6.1 Limits from this analysis 

Here we present the study of the anomalous magnetic and electric dipole moments of 
the tau lepton based on the analysis of the e+e~ -^ e^e~T~^T~ cross-section. The study 
of anomalous couplings of tau leptons to photons at LEP in this channel was proposed 
in [18]. 

To model the contribution of non-SM anomalous magnetic and dipole moments we 
use the calculation by Cornet and Diana [19]. The calculation is based on computation 
of the matrix element of the process 77 -^ t^t^ in leading order of QED and its trans- 
lation to the cross-section of the e"'"e~ — > e"'"e~r+r^ process using the Equivalent Photon 
Approximation (EPA) [20]. The EPA parameter {—q^)max (the upper limit of the integra- 
tion over 4-momenta of the emitted photon) was chosen such that the total cross-section 
predicted by EPA (with SM values of anomalous electromagnetic moments) agreed with 
BDKRC calculation. According to the calculations [19] each of the anomalous terms 
of (5) would mainly modify the rate of tau pair production in the barrel region of the 
detector where the experimental selection has largest efficiency. This leads to a larger 
selection efficiency for the anomalous term contribution, improving in principle the limits 
obtained on anomalous moments. However in this paper we conservatively assumed that 
the standard and anomalous contributions have the same selection efficiency. 

Figure 7 shows how the total cross-section changes as a function of the anomalous 
magnetic moment and as a function of the electric dipole moment. The three lines on 
each plot represent the calculation with ^/s =182.7, 195.5 and 205.0 GeV. Increasing 
the collision energy slowly increases both non-SM contributions. However, increasing 
the magnitude of the anomalous magnetic moment can either increase or decrease the 
cross-section while increasing that of the electric dipole moment tends only to increase 
the cross-section. 

To compare the experimentally measured values of the cross-sections to the non-SM 
calculation, they were first converted from accepted to the total cross-sections, taking into 
account the 1.44% difference due to the signal definition (see section 2). The validity of 
applying SM conversion factors is supported by the fact that the measured cross-sections 
are in good agreement with the SM prediction, which guarantees the smallness of the 
non-SM contribution, and by the fact that the correction itself is small. 

Fits to the cross-sections measured in 1997, 1998, 1999 and 2000 were performed taking 
ttr and dj. as parameters. When fitting for a,-, the value of dr was set to its SM value and 
vice versa. The errors on the cross section measurements were taken as the statistical 
and systematic errors added in quadrature. 

To quote the obtained limits we used the following convention: 

/L r-co -j^ _ Qj^ 

exp (-xV2) dar = / exp (-x'/S) ^Or = ^7^ (8) 

-00 J R ^ 

where CL is the desired confidence level and L and R are lower and upper limits. A 
similar definition was used for d^-- We quote central values /i and errors a for moments 
according to 

where R and L are calculated with 68.3% confidence level. 

Figure 8 shows the x^ ^^ ^ function of the anomalous magnetic moment and as a 
function of the electric dipole moment. The results of the fit are: 

-0.052 <ar< 0.013, 95% CL 

\dr\ < 3.7 • 10-16 e . cm, 95% CL. 



11 



The limit on a^. improves the current PDG hmit [21] based on the L3 result [14]. 

Figure 9 shows the the measured cross- sect ion, average LEP2 cross-section and SM 
expectation as a function of ^/s. Two bands superimposed on the plot represent the 
allowed region for the cross-section variation due to anomalous magnetic and electric 
dipole moments. The results expressed in the form of central value and error are the 
following: 



ar = -0.018 ±0.017, 
dr = (0.0 ±2.0) -10-16 e-cm. 



7 Conclusion 

We have studied the reaction e"'"e- -^ e^e~T^T^ with the data collected with the 
DELPHI detector during LEP2 operation in the years 1997-2000. The average LEP2 
cross-section was found to be 429±17 pb compared to 447.7 pb expected from the Stan- 
dard model. The measured/predicted ratio 0.96±0.04 agrees with the QED prediction 
at the level of one standard deviation. The measured cross-sections were used to extract 
limits on the anomalous magnetic and electric dipole moments of the tau lepton. The 
95% CL limits obtained are 

-0.052 <ar< 0.013 
\dr\ < 3.7- 10-16 e-cm. 
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Figure 2: The distributions of invariant mass, total energy, transverse energy and trans- 
verse momentum of the pair of charged particles. Preselection cuts are applied. The 
points are 1999 data, the open histogram is the simulation of background processes and 
the shaded histogram is the simulation of e"'"e~ — >■ e~^e~T~^T~ events. The simulation is 
not corrected for the trigger efficiency. 
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Figure 3: Distribution of the dE/dx pull for the electron hypothesis (left) and for the 
muon hypothesis (right) with all selection cuts applied except the cut on the variable 
shown. Points are 1999 data, the open histogram is the background, and the shaded 
histogram shows the e"'"e~ — > e"'"e~r"'"r~ signal events. The simulation is not corrected 
for the trigger efficiency. 
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Figure 4: Left: specific energy loss as a function of particle momentum for electron candi- 
dates after preselection cuts and the cuts on the electron and muon pulls. The horizontal 
line shows the first cut against kaons and protons. The points are 1999 data. Right: 
distribution of dE/dx pull for proton hypothesis for electron candidates. The hatched 
histogram is the background from e"'"e~ — >■ e"'"e~qq events and the shaded histogram is 
the rest of the background. The cuts against protons are indicated by arrows. All other 
selection cuts are applied. The points are 1999 data. 
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Figure 5: Momentum distribution for electron candidates (left) and non-electron can- 
didates (right) for selected events from 1997-2000 data. The distributions of simulated 
events are corrected for trigger efficiency. The dip in the electron momentum distribution 
near 1 GeV/c is caused by the cut against protons: the electron and proton dE/dx are 
equal in this region of momentum. 
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Figure 6: Visible invariant mass distribution for selected events for combined 1997-2000 
data. The distribution of simulated events is corrected for trigger efficiency. The mass 
was calculated using all detected charged particles and photons. The simulation was 
corrected for the trigger efficiency. 
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Figure 7: Total cross-section change as a function of anomalous magnetic moment and 
as a function of electric dipole moment. 
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Figure 8: x^ ^^ ^ function of anomalous magnetic moment and as a function of electric 
dipole moment. 



17 



-Si 



H 
H 

^ 
^ 






600 



550 



500 



450 



400 



350 



30t 



DELPHI 

• Measurements 

■ Mean 

- Standard Model 




a allowed (95% CL) 

d^ allowed (95% CL) 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



70 175 180 185 190 195 200 205 210 215 220 

^s (GeV) 



Figure 9: Measured cross-section (circles), average LEP2 cross-section (square) and SM 
expectation as a function of ^/s. The two bands show the cross-section variation allowed 
due to anomalous magnetic and electric dipole moments within 95% limits from this 
analysis. 



